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Abstract. The relativistic properties of the three-nucleon system are investigated using 
the Eaddeev equations within the Bethe-Salpeter approach. The nucleon-nucleon inter¬ 
action is chosen in a separable form. The Gauss quadrature method is used to calculate 
the integrals. The system of the integral equations are solving by iterations method. The 
binding energy and the partial-wave amplitudes (*5o and ^5i) of the triton are found. 


1 Introduction 

Three-body calculations in nuclear physics are of great interest used describing three-nucleon bound 
states (^He, T), processes of elastic, inelastic and deep inelastic scattering the leptons off light nuclei 
and also the hadron-deuteron reactions (for example, pd —» pd, pd —» ppn). Study of the nuclei 
^He and T is also interesting because it allows us to investigate further (in addition to the case of the 
deuteron) evolution of the bound nucleon thereby contributing to the explanation of so-called EMC- 
effect. In quantum mechanics the Eaddeev equations are commonly used to describe the three-particle 
systems. The main feature of Eaddeev equations is that all particles interact through a pair potential. 

However at the high momentum transfer relativistic effects should be taken into account. The 
Bethe-Salpeter (BS) ||T] equation is one of the most consistent approaches to describe the NN inter¬ 
action. In this formalism, one has to deal with a system of nontrivial integral equations for both the 
NN scattered states and the bound state - the deuteron. To solve a system of integral equations, it is 
convenient to use a separable Ansatz ll for the interaction kernel in the BS equation. In this case, 
one can transform integral equations into a system of algebraic linear ones which is easy to solve. 
Parameters of the interaction kernel are found from an analysis of the phase shifts and inelasticity, 
low-energy parameters and deuteron properties (binding energy, moments, etc.). 

The relativistic three-particle systems are described by the Eaddeev equations within the BS ap¬ 
proach - so called Bethe-Salpeter-Eaddeev equations. In the paper the three-nucleon nuclei is studied 
with the simplest Yamaguchi separable potential lO. All nucleons have equal masses and the scalar 
propagators instead of spinor ones are used for simplicity. The spin-isospin structure of the nucleons 
is taken into account by using the so-called recoupling-coefficient matrix. The work mainly follows 
the ideas of the article a. 
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The paper is organized as following: in Sec. 2 the two-particle problem is considered, in Sec. 3 - 
three-particle equations. In Sec. 4 the calculations and results are given. The summary is in the Sec. 5. 


2 Two particle case 

Since the formalism of the Faddeev equations is based on the properties of the pair nucleon-nucleon 
interaction only some conclusions of two-body problem are given here. 

The Bethe-Salpeter equation for the relativistic two-particle system is taken in the following form: 

T{p,p'-,s)^V{p,p')+ Jd‘^kV(p,k)G(k-,s)T(k,p';s) (1) 

where T{p,p'\ s) is the two-particle T matrix and V(p,p') - kernel (potential) of the nucleon-nucleon 
interaction. The free two-particle Green function G(k; s) is expressed, for simplicity, thought the 
scalar propagator of the nucleons 

G-\k-, s) = [(P/2 + kf - ml + ie][(P/2 - kf - ml + ie\. (2) 

To solve equation Q the separable Ansafz for the nucleon-nucleon potential V(p, p') is used (rank- 
one) 

V(po,p, Po, p') = Agipo, p)g(p'^, p'). (3) 

In this case the two-particle T matrix has the following simple form: 

T(po,p,Po,p'-,s) = T(s)g(po, p)g(pQ, p') (4) 


where 


[t(^)]-' = T-' 


• ^OO ^OC 


k^dkg^(k°,k)G(k°,k; s). 


As a simplest assumption the relativistic Yamaguchi-typ& form factor gY(po,p) is used 

1 

griPa, P) 


-Po + P^+0- 


:2’ 


( 5 ) 

(6) 


with parameters (3 and A are chosen to describe experimental data (deuteron binding energy, low- 
energy scattering parameter and phase shifts). The values of the parameters are given in Table[T] 

Table 1. Parameters for T(i) and gy(po, p) for S-states 


Parameter 

‘5o 


A (GeV'*) 

-1.12087 

-3.15480 

/3 (GeV) 

0.287614 

0.279731 


3 Three-particle case 


The relativistic three-particle system can be described by the Bethe-Salpeter-Faddeev equations 
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where fall matrix T - Gi is the two-particle (j and n) Green function (/jn is cyclic permu¬ 

tation of (1,2,3)): 


Gi(kj, ^„) = 1 /{k^j -m\ + ie)l(kl - ml + ie), ( 8 ) 

and Ti is the two-particle T matrix which can be written as following: 

TAkuki, h-, k[,k'^, k'^) = - K'.)Ukj, k„-k], k'J. (9) 

For the system of equal-mass particles the Jacobi momenta can be written in the following form: 

( 10 ) 

Using expressions (10 1 the equation Q can be rewritten as 


Pi = ~ ^n), qi ki, K^ki+k2+ k^. 

2 3 


T^‘\pi, qr, p'i, ql s) = { 27 T)'^ 6 ^'^\qi - q'i)Ti{pi; p'; s) 
dp'i 


( 11 ) 


' f k'l)[T^J\p'J, q''; p\, qls) + T^^Xp”, q”\ p', q'f, i)] 


It is suitable to introduce the amplitude 'I'^'l(p;, q',; s) for the bound state of the three particles in 
the following form: 


^^‘\pi,qr, s) = <p,-,^,|r®|MB) = ^'LM{p,q\ s). 


( 12 ) 


where Mg - ^ - 3mf^ - Eg is the mass of the bound state (triton) and s - is the total momentum 

squared. For the equal-mass case all functions are equal to each other. 

The total orbital angular momenta of the triton can be presented as L = l+A, where / is the angular 
momentum corresponding to nucleon pair with relative impulse p and A is the angular momentum 
corresponding to relative impulse q. 

To separate the angular dependence the amplitude can be written in the following form: 

'i'LMip, q\s)^Yj IPl’ ^0’ I'll; ")‘^Sm(P- q)> (13) 

aA 


with the angular part 

mfi 


(14) 


where the two-nucleon state with spin s, angular moment / and total momentum j (a = is 

introduced, C is the Clebsch-Gordan coefficients and Y is spherical harmonics. 

If one consider the rank-one separable two-nucleon interaction the amplitude can be written 
as 


'*' 2 (^ 0 , s) = g''‘‘\pQ,py‘^’(sW°l(qQ,q-, s). 


s(a)/ 


where function satisfies the following integral equation: 




• _ ^00 r^QC 

yx Zdq'ZZ? 


'}l^\qo,q; qQ,q';s) 

r(“')[(fV^ + ^') 2 -^' 2 ] 

— %/(qo,q;s). 


(5 a/s - q'o)^ - q^ -mA + ie 


(15) 


(16) 
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with 


q-, q'o, q'-, s) = C^aa') J d cos §qq,Kf^'l(q’ *^ 99 ') 
g^^K-^qo - q'o, lq/2 + q'l) 0 ^‘''*(q'o + ^'q, Iq + q72|) 


(5 aAs + g'o + q'o)^ - (q + q')^ - 


+ le 


and 


Kffliq, q, cos d-qq') = (47t) 


3/2 V2/1 + 1 


2L+ 1 

I] CfZoC^Z'A'm-nr'Y^imid, 0)Yr„-(d', 0), 0) 

mm' 


(17) 


(18) 


where ^ ^ 

cos?? = (| +q''cos?7,^0/l| +q'l, cost?' ^ (q + |-cos i?^,0/lq + yl 

and C(aa’) is the spin-isospin recoupling-coefficient matrix. 

Considering the ground state of the three-particles state L - 0 and assuming the internal orbital 
angular momenta are equal to zero I - A - 0 one should take into into account only two states: 
a = (*5o,^ S i). In this case function - 1 and system of integral equations reads 


O^“i(q' 0 , ^ ^ dgoj q'^dq Z‘'“°'\qQ, q, q'^, q\ s) 

+ q'^f - q'^] 


(19) 


(\^rs-q'oY-q'^-^ 


-Q>^'^'\q'^,q'-s). 


le 


Here the so-called effective energy-dependent potential Z is 
Z^‘“‘’\qo, q, q'o, q'; s) = C^aa') J" d{cosdqq,) 

g^‘'^i-\qo -q'Q,\-\fi -+ \q'Q ,Iq + ^q'l) 


( 20 ) 


(5 Vi + q'o + g'o)^ - (q + q')^ 


with 


Q 


{aa') 


1 

3 

4 


_3 

V 

4 


( 21 ) 


The system of integral equations ( [T9p0 1 has the number of singularities, however in the case of 
the bound three-particle system { yfs < hniN) all this singularities do not cross the path of integration 
on qo and thus do not affect to the Wick-rotation procedure qo ^ iq 4 . 

The system of Eqs. ( [T9p0| l after the Wick-rotation procedure is well analytically defined and can 
be solved by different numerical methods. One of them is discussed in the next section. 


4 Solution and results 

In order to solve the system of integral equations the iterations method is used. The mappings for the 
variable of integration on q [0, 00 ) and q 4 (- 00 , cm) to [-1,1] interval are introduced. 
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The binding energy of the three-nucleon system is satisfied the following condition (see details 
in 13): 


lim 

rt—>oo 


o«(5) I 


= 1 


( 22 ) 


where n is number of iteration. 

Result of calculations for the binding energy is Eb = 10.99 MeV which should be compared to the 
experimental value 8.48 MeV. The difference can be explained by the simplicity of the using separable 
kernel (rank-one) of the nucleon-nucleon interaction. 

The obtained partial-wave amplitudes are shown in the Figs. m- 


Re 'Sij-function Im *Su-function 




Figure 1. Real (left) and imaginary (right) parts of the *So partial-wave state as a functions of at different q 
values: solid black line - for q = Q Fm“*, dashed red line - for q = 0.5 Fm“' and dotted-dashed blue line - for 
^ = 1 Fm ' 


Re ^Sj-funetion Im ^Sj-funetion 




Figure 2. Real (left) and imaginary (right) parts of the ^5 1 partial-wave state as a functions of q^ at different q 
values: solid black line - for q = Q Fm“', dashed red line - for q = 0.5 Fm“’ and dotted-dashed blue line - for 
q = \ Fm ' 
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Re *Sj-function Re ^Sj-funetion 




Figure 3. Real parts of the *5o(left) and ^5 1 (right) partial-wave state as a functions of q at different q^ values: 
solid black line - for ^4 = 0 Fm“', dashed red line - for ^4 = 0.1 Fm“* and dotted-dashed blue line - for q^ = 0.2 
Fm-' 


5 Summary 

In the paper three-body system is investigated by using Bethe-Salpeter-Faddeev equations. The 
nucleon-nucleon interaction is taken in the rank-one separable form. The relativistic generalization 
of Yamaguchi-type functions are chosen in the calculations. The parameters of the nucleon-nucleon 
potential for the '5o and \ partial-wave states reproduce the low-energy scattering parameter and 
deuteron properties as well as phase shifts up to the laboratory energy 100-120 MeV. The BSF integral 
equations are solved by using iterations method. The binding energy of the triton and amplitudes of 
the '5o and i partial-wave states of the triton are found. 

The large overestimation of the binding energy of the triton is found. To improve results the rank 
of the separable kernel should be increased. Also another partial-wave states, such as P- and D-states, 
and spinor propagators for nucleons should be taken into account. 


Acknowledgments 

One of the authors (S.G. Bondarenko) thanks organizers of the International Conference “Mathemat¬ 
ical Modeling and Computational Physics, 2015” (Stara Lesna, Slovakia, July 13 — 17, 2015) and 
personally Professor Michal Hnatic for invitation, support and hospitality. 


References 

[1] E.E. Salpeter, H.A. Bethe, Phys. Rev. 84, 1232 (1951). 

[2] S.G. Bondarenko, V.V. Burov, A.V. Molochkov, G.I. Smirnov, H. Toki, Prog. Part. Nucl. Phys. 
48, 449 (2002). 

[3] Y. Yamaguchi, Phys. Rev. 95, 1628 (1954); Y. Yamaguchi, Y. Yamaguchi, Phys. Rev. 95, 1635 
(1954). 

[4] G. Rupp and J. A. Tjon, Phys. Rev. C 37, 1729 (1988). 

[5] R. A. Malfliet and J. A. Tjon, Nucl. Phys. A127,161-168 (1969) 

















